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ABSTRACT 


The radio images reconstructed by using Synthetic 
Aperture Radio Telescopes are distorted due to the effect of 
ionospheric Irregularities. The effects are especially pronounced 
at lower frequencies. The effects have been studied for the Giant 
Metre Wave Radio Telescope (GMRT) being built at Pune, India. 
GMRT being a metre wave radio telescope Is very susceptible to 
these Image distortions due to the effect of the ionospheric 
Irregularities. A new approach to integrate the work of both the 
astronomers and the propagation community has been tried using 
the mutual coherence function (MCF). Simulations to study the 
effects of different parameters like frequency, declination, sun 
spot number etc. have been carried out. The Ionospheric 
Irregularities have been modeled by the UBMOD program for 
calculation of scintillation Indices developed by Fremouw. The 
effects can be clearly seen In the color coded photographs 


attached . 
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CHAPTER - 1 


INTRODUCTION 

Astronomy is traditionally "looking at the sky" and 
interpreting the message that starlight brings. But objects in 
space emit more than ordinary light. They radiate other wave- 
lengths (X-rays, ultraviolet, infrared and radio waves) and most 
of the exciting astronomical discoveries of the past 20 years 
quasars, pulsars, black holes, proof of the Big Bang have been 
made by studying these other radiations. 

1.1 RADIO ASTRONOMY 

Radio Astronomy is the "Grand old man” of the new 
astronomies. X-ray, ultraviolet and infrared radiations are 
absorbed by Earth’s atmosphere, and it is difficult to study them 
except from high mountains, balloons, rockets or preferably 
satellites completely above the atmosphere. Radio waves from 
space do reach the ground, however, so their observations is an 
easier proposition. Nearly 60 years have passed since cosmic 
radio waves were first detected, and in that time radio 
astronomers have overtaken their optical colleagues in many ways. 


Radio telescopes can 

now 

"see" details 

a thousand 

times finer 

than optical astronomers 

can photograph 

and 

the 

positions 

of 

radio sources can be 

pin 

pointed to five 

times 

the 

accuracy 

of 


optical positions in the sky. 
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The two major criteria o£ performance of a telescope are its 
sensitivity and its resolving power. The first function 
determines its ability to observe objects otherwise too faint to 
be seen with the unaided eye. The second function determines 
telescope’s ability to separate sources that would appear as a 
single source to the naked eye. The theoretical resolving power 
of a telescope, defined as the smallest angular separation 
between two sources which can be seen as two distinct sources, is 
proportional to X/D, where X is the wavelength of the radiation 
and D is diameter of the aperture. 

To obtain a resolution equivalent to that of a 200 inch 
optical telescope, the diameter of the collecting surface for a 
radio telescope would be some 150 miles at 1000 MHz. 

The skillful applications of the techniques of radio inter- 
ferometry coupled with the evolution of computer technology, has 
led to resolving powers that were inconceivable in the early days 
of radio astronomy. In the past few years radio astronomers have 
obtained resolutions of 0.0001 arc seconds which constitutes an 
improvement of some four orders of magnitude over the resolution 
of earth based optical telescopes. 

1.2 APERTURE SYNTHESIS AND GMRT 

The high resolutions of 0.0001 arc second does not violate 
the relation X/D. The requirement of making the aperture D very 
large is accomplished by "aperture synthesis”, a technique in 
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which the radio waves collected by two (or more) instruments upto 
tens, hundreds and even thousands of kilometers apart are 
recorded simultaneously and subsequently added together. The 
addition must be done, of course, in a special way; it is 
necessary to arrange the observations precisely so that the phase 
relations of radio waves arriving at the telescope participating 
in the synthesis is preserved. 

Since the earth rotates, the relative position of two radio 
telescopes is continuously changing. As a result the signal 
intensity at the output of the Interferometer passes rapidly 
through a succession of maxima and minima called interference 
fringes . 

The actual shape of the intensity variations is a sine wave. 
The amplitude and phase of the sine wave embody information about 
the structure of the celestial object on a certain angular scale, 
which depends on the distance between two telescopes as seen from 
the source. If the distance between the telescopes is changed, 
one gets a different pattern, which then contains information 
about the structure of the object on a different angular scale. 

It Is therefore possible to measure the fringe pattern for a 
range of telescope seperations and so to cover a range of angular 
scales on the sky. The amplitudes and phases of the fringes, 
sampled for a large number of seperation of the two telescopes 
are called the visibility function. 
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This visibility function happens to be the Fourier Transform 
of the brightness distribution of the source in the sky. Hence if 
the visibility function is measured for all possible baselines, 
then, theoretically at least, the brightness distribution can be 
obtained by Fourier transformation. 

In practice it Is possible to obtain the fringe pattern for 
all seperat ion of two telescopes out to a maximum seperat ion , D. 
This enables one to obtain the same image that one would get from 
a single giant telescope of diameter D. The different baseline 
separations D can be obtained either by moving one of the 
antennas or by having large number of different antennas. Uhen 
the rotation of the earth is exploited to increase the number of 
separations, the resulting operation is called earth-rotation 
aperture synthesis or super synthesis. Most of the large radio 
telescopes built in the last two decades follow this principle. 
In India, the Giant Metre-Uave Radio Telescope (GMRT) is being 
setup at Pune. GMRT will consist of 30 fully steerable 45 m 
diameter parabolic dishes and is being designed to operate in six 
separate frequency bands between 38 MHz and 1420 MHz. The GMRT 
will operate as an "Earth Rotation Synthesis Radio Telescope". 
Its 30 antennas will form 435 interferometric pairs with the 
length and orientation of the projected baselines changing 
continuously with the diurnal rotation of the earth. Uhat makes 
GMRT unique is the relatively lower frequency range of operation. 
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1.3 PROPAGATION EFFECTS 

The fundamental measurements in the aperture synthesis 
telescope are the amplitude and phase of the fringe visibility. 
Both of these parameters are affected by the passage of the wave 
through the earth’s atmosphere. The irregular structure of 
refractive index of the atmosphere causes random variation in 
amplitude and phase called scintillation. At lower frequencies 
the scintillation due to the ionosphere are the predominant 
factors. Since Pune is In the equatorial F-region ionospheric 
scintillation region, the study of ionospheric scintillation on 
the reconstructed images becomes Important. This thesis tries to 
address this problem. 

Uhen a wave propagates through the Ionosphere, there are two 
major effects that come into play 

I) The regular effects, and 

II) The irregular effects. 

The regular effects are primarily 

a) Refraction, due to change in permittivity which changes 
regularly with the electron density in the different regions 
of the ionosphere 

b) Absorption 

c) Faraday rotation which changes the plane of polarisation of 
the radio wave, and 

d) Introduction of an excess path length. 
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All these changes are related to the change In permittivity 
of the ionosphere given by 

N N ■* electron density 

. e 6 

&£ <* — 

u v ■* frequency 

-2 

Clearly all these effects scale as v and their effects are 
well studied. Davies (1990) has given an excellent treatment of 
ionospheric propagation and reviews of particular relevance are 
given in (Evans, 1968) and (Hagfors, 1976). 

The irregular effects a 

The irregularities in the ionosphere perturb the cosmic 
signal. The predominant scale sizes in the ionization 
irregularities are found to be a few kilometres or less. 

The effects of these fluctuations have been studied 
extensively at frequencies between 20 and 200 MHz. Hewlsh (1952), 
Booker (1958) and Lawrence, Little and Chivers (1964) have 
reviewed the early results and techniques. A comprehensive review 
of theory and observations of ionospheric fluctuations can be 
found in (Crane, 1977), (Fejer and Kelley, 1980) and (Yeh and 
Liu, 1982) and a summary of global morphology can be found in 
(Aarons, 1982). 

The most important effect of ionospheric irregularities, is 
a random fluctuation in phase and amplitude known as phase 
scintillation and amplitude scintillation respectively. Spoelstra 
& Kelder (1982) have studied the effects of ionospheric 
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scintillation on synthesis telescopes. Due to change in phase of 
the measured visibility function, the image can be distorted as 
the phase carries the information about the spacing of the point 
sources in the image. 

1.4 SUMMARY OF PRESENT WORK 

In this thesis a new approach to study the effects of 
ionospheric scintillation on the reconstructed image has been 
tried, with particular emphasis on GMRT . The present work tries 
to integrate the studies made by astronomers and the ionospheric 
propagation community. The common element and the starting point 
is the Mutual Coherence function or the Spatial Coherence 
function of Fomalont (1978). A lot of studies have been carried 
out on the effect of ionosphere on the mutual coherence function 
with the most recent and relevant one being that done by 
Bhattacharya & Rastogi (1991). Although they have carried out 
studies for single telescopes, the results are of a more 
fundamental nature and have been extended to spaced receivers, 
working as interferometric pairs. 

The thesis is organised into five chapters. The second 
chapter describes the process of image reconstruction from 
visibility measurements and chapter 3 gives the details of 
ionospheric scintillation and its effect on the Mutual Coherence 
function. Chapter 4 discusses the details of the computer code 
generated for the work and chapter 5 discusses the results. 



CHAPTER - 2 


INTERFEROMETRY AND RADIO IMAGE MAPPING 


2. 1 INTERFEROMETRY THEORY 


The theory o£ aperture syntheslc has been well summarised by 
Thomson, Moran and Swenson (1986). A simple interf erometeric 
arrangement is shown in Fig. 2.1. 


For a point source, each antenna delivers the same signal 
voltage V with a phase difference determined by the baseline D 
and the source direction O. The integrator has a time constant 
2T. The output of the integrator resulting from a point source is 


r(r ) 


1 

2T 


T 

J V(t) V(t - t) dt 
-T 


( 2 . 1 ) 


The integration time 2T is chosen to be very large as compared to 
Av \ where Av is the finite bandwidth of the two pre-detection 
amplifiers. Clearly the eqn. (2.1) is an autocorrelation 
function. 


The time delay between the outputs of the two antennas is 

given by r = (D/C) sin#, r is called the geometrical delay and 
g & 

C is the velocity of light. If V(t) has a sinusoidal variation of 
the form V(t) = sin2rci>t. The output of the multiplier then is 


proportional to 
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F = 2sin(2m>t) sin2nv(t-T ) = cos2m>r - cos(4rci>t) cos(2»zvt ) 

fi fi fi 

- sin(4ni>t) sin(2rrvr ) . 

£ 

By filtering the more rapidly varying terms in F can be 
filtered out leaving 

f 2nD . . 1 f 2nDZ T 

F = COS I BtrtS I = COS I — > I 

where £ = sin# is the direction cosine measured with respect to 
the baseline. 

The power spectrum of a signal is the Fourier transform of 
the autocorrelation function of the signal. This is known as the 
Ui ener-Khinchin relation. This can be written mathematicaly as 

oo 

r(r) e ~ j2nVT dr = |H(i>) 2 | (2.2a) 

-to 

and 

oo 

|H(v) 2 | &l> = t(t) (2.2b) 

-oo 

2 

where H(i>) is the amplitude (voltage) response, and hence |H(v)| 
is the power spectrum of the signal input to the multiplier. If 
we consider the cosmic signal to have a constant amplitude 
spectrum (which is valid almost universally), the amplitude 
spectrum H(v) is determined solely by the passband 

characteristics (frequency response) of the amplifiers. Thus the 
output of the interferometer as a function of the time delay t is 
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the Fourier transform of the power spectrum of the cosmic signal 
as bandlimited by the amplifiers. 

If V( t ) has a sinusoidal variation of the form sin(2nyt) , 
then the output of the multiplier and the Integrator comes out to 
be of the form 


r(r) = cos 


2n D? y 1 
o 


X Bandwidth pattern. 


where % = SinO is the direction cosine measured w.r.t. the 

baseline 

= Center frequency of the band limited amplifier, and 

2 

Bandwidth pattern depends upon the form of |H(v)| . For a uniform 

power spectral density, the bandwidth pattern = -r^— - so 


nDfAv/C 


that the correlation response is given by 


2Tlv D 


, . f o . . 1 sin (n D 

r(T) = =»[ — — sine j 


K A i>/C) 


£ Av>/C 


Figure 2.2 shows the point source response of the type given 
by the above equation. 

In mapping applications, observations are made in the 

vicinity of the bandwidth pattern maxima and this condition is 

achieved by changing the instrumental delay r. periodically so as 

to keep t = r - t. suitably small, 
g i 

The instrumental delay is calculated for a nominal source 
position & which becomes the reference position. Consider the 
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response to a narrow source at position P in the vicinity of 

Let 6 - & - &* where S’ is very small. Then, interferometer 
o 

response is given by 

F(£ ’ ) = cos ( 2n u £ ') 

f V o D 1 

where £’ = sin<9 at &' is the direction cosine, u = | — | cob& q = 

^D/X q J co8^ o is a spatial frequency term which is measured in 
cycles per radian, since the spatial variable £ ’ being small, can 
be measured in radians. 

2. 2 SOURCE SYNTHESIS 

If B^(£’) is the brightness of the source, we obtain the 
interferometric output as a function of £ ’ . 

R(£ ’ ) = Ai> J cos J^2nu (£’-£" )J A(£"-£ ’ ) B^”) d£” 

source 

where A(£ ’ ) is the antenna response variation with direction. 

Now for a particular value u = the Fourier transform of 

the fringe term cos(2nu o £ ’ ) is given by 

cos^2nu o £’J = ~ J^SCu-u^) + «5(u+u o )J 

Defining 

bj^Cu) 4» B x (£ ’ ) and r(u) 4» R(£ ) 

and assuming broad enough antenna beams so that gain is constant 
over the source we have 
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Fig. 2.3 : Power reception pattern of a typical antenna, which is 
a polar plot of effective area versus angle of arrival 
of plane-wave radiation. 
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r(u) = j AAi> j^b 1 (u o ) <5 (u-u q ) + b^-u^ «5 (u+u q )J 

Hence r(u) constitutes the angular (spatial) frequency spectrum 

of the source brightness function as filtered by the 

interferometer, which passes only the spatial frequencies ± u . 

o 

Consider an interferometer with an east-west baseline and an 
extended cosmic source that passes thrugh the interferometer 
reception pattern from east to west as the earth rotates. Each 
hour angle of the source corresponds to a different value of the 
angle between the baseline and the direction to the source, and 
thus to a different spatial frequency u. The variable also 

increases with time, so the output of the interferometer varies 
cosinusoidal ly (or sinusoidally) with time over small intervals. 
The earlier equation shows that at any instant the amplitude of 
the fringe pattern from the correlator r(£) is equal to the even 
part of the Fourier transform of B (£ ' ) , that is, to the cosine 
transform of the brightness with respect to the point on the 
source at position angle at that time. Uith respect to this 
same point on the source at an instant one-quarter of a fringe 
cycle later in time the output fringe pattern is equal to the 
sine transform, which is the imaginary part of the brightness 
Fourier transform. Thus the output fringe pattern contains both 
the real and imaginary parts of b^fu), the Fourier transform 
(£ ’ ) . In practice, a phase reference position on the sky is 
designated, central to the region of the Interest, and the fringe 
pattern is then computed for a hypothetical point source of unit 
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flux density at this position. The amplitude and phase of the 
observed fringe pattern are measured with respect to the 

hypothetical fringes. The resulting values, when expressed as a 
complex quantity, are known as the complex visibility V(u), which 
is equal to b^(u). The complex visibility, as derived from the 
measured output of an interferometer, includes the effects of the 
antenna beams, the finite receiving bandwidth and any 

instrumental imperfections in addition to the effects of the 
atmosphere. In practice the antennas often track the source, so 
is multiplied by A^ . Then 

R(?’) = J cos [2nu (K ’ ~ ?")J A «f") dZ” 

source 

If the phases and amplitudes of the fringes are observed at 
frequent enough intervals in u, the complex spatial spectrum 
represented by V ( u ) is obtained. Then by Fourier transforming 
this spectrum and correcting for A^ the source brightness 
distribution can be reconstructed. 

Uhl 1 e observations of a source for many hours with a single 

baseline can yield, in principle, a wide range of spatial 

frequencies, in practice they usually cannot produce a complete 

enough spectrum to permit satisfactory reconstruction of a source 

profile. Several different baseline are thus required. These may 

be obtained from n (n -l)/2 pairs in a set of n 

a a u 

form an array. 


antennas that 
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2.3 TWO DIMENSIONAL MAPPING 


In the particular case of an east-west baseline and an 
equatorial source, the source can be mapped in the east-west 
dimension only, as the spatial frequency spetrum so derived has 
only one frequency dimension, u. To synthesize a map in two 
dimensions requires a two dimensional spatial frequency spectrum 
in the (u,v) frequency domain, v being the north-south frequency 
component. This may be accomplished for an equatorial source by 
observing with a two-dimensional array of interferometers. The 
-two-dimensional situation is similar to the one— dimensional one 
but more complicated. It can be understood by using the more 
general form of the visibility function which is given as 

j = J A^o) B(>) e ” dO . 

4jt 


where V 


A n O) 

B<>) 

dft 


complex fringe visibility 
baseline vector of the two antennas 
is the normalized antenna area 
is the radio brightness 

an element of solid angle on the sky & — S — 
ahown in Fig. 2.4. 


(.K , D, O are the direction cosines of s. S q is the reference 
direction and its direction cosines are (0,0,1). 


D X .S » U* + V7) + wC 


Since 
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For small source 


Since 


<y = S - S is small 
o 


£ , n are small and £ — 1 • 


Under these conditions we have 


dCi = * d y =* dZ dn 


A 72 2 

/ 1 - - T) 


and we have 


CO 00 f 

u,v,w) = J | A N (?, D) B(? , Y) ) x exp | - j 2n ju{ 


V( 


+ V 7 ) 


-00 -00 


+ W(C - 1) 


ll dy > 

J J /^v 


which reduces to 


oo oo 


V(u, v,w) = V(u , v , 0 ) = J I A n (?,yj) B(?,r>) e 3 


- j2fi(u£ + vt >) 


dZ dr, 


-oo -co 


2.4, ANTENNA SPACING COORDINATES AND C u, v'J LOCH 

A right handed Cartesian coordinate system in which X and Y 

are measured in plane parallel to the earth’s equator, X in the 

meridian plane (defined as the plane through the poles of the 

earth and the reference point in the array), Y towards the east, 

and Z is measured toward the north pole. If (X^, are t * ie 

components of 5^ in the (X, Y, Z) system, the components (u,v,w) 
A. 


are given by 



Spacing vector locus in the (u # v) plane. Tne 
lower arc in the diagram represents the locus 
of conjugate values of visibility. Unless the 
source is circumpolar the cutoff at the 
horizon limits the lengths of the arcs. horizon 
limits the lengths of the arcs 
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Fig. 2.6 : The coordinate system for the baseline measurement 
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u 


sinH 

cos H 

0 


X. 







X 

v 

= 

-sin<5 cosH 

sin<5 sinH 

cos<5 


Y x 

w 


cos6 cosH 

-cos6 sinH 

sin6 



- 


- 






where (H,<5) are the hour angle and the declination of the phase 
reference position. Examination of above equations shows that 
v 


u 2 + 


Z. cos<5 "1 
X o 

sin<5 

° J 


+ < 


which defines an ellipse in the (u,v) plane with hour angle as 

the variable. The ellipse is centred on the v axis at (u,v) = 

2 2 

X x+ Y x 

/ “5 2 * 

X. + Y., . Since V(-u, -v) = V (u,v), any observation 
XX 

supplies simultaneous measurements on two arcs, which are part of 
the same ellipse only if Z^Q . Fig. 2.6 describes the co-ordinate 
system used. 


2.5 DISCRETE FOURIER TRANSFORMATION 

The discrete transform is very widely used in synthesis 
mapping because of computational advantages. With the discrete 
transform the function V(u,v) and B(£ , y?) are expressed as 
rectangular matrices of sampled values at uniform Increments in 
the two variables involved. The main advantage of this system is 
that the Fast Fourier transform ( FFT ) algorithm can be used to 
reduce the computation time. The two-dimensional form of the 
discrete transform for a fourier pair f and g is defined by 
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M-l N-l 

_ 1 V V „rt, , N _ j 2nrkm/M -j2nln/N 
f(m,n) ~ fjfj > ) e(k,l) e e 

k=0 t^O 

and the inverse is 

«(k.l) = I I £ (m , n ) e ^ ^rckm/M e j 2rc 1 n /N 

m=0 n=0 


To apply the discrete transform to the synthesis-mapping 
problem we first use an interpolation to obtain values of V(u,v) 
at points separated by Au in u and Av in v as shown in Fig. 2.7. 

The dimensions of the (u,v) plane that contain these data 
are MAu by NAv . In the (£ , y) ) plane the points are spaced A£ in £ 
and Ay) in y) and the map dimensions are MA£ by NAy) . The dimensions 
in the two domains are related by 


Au = 



-1 


Av 



-1 


A K = 

Ay) = 



The transformation between V(u,v) and B(£ ,T )) is obtained by 
substituting g(k , i ) = B(kA£ , lAy )) and f(m,n) = V(mAu f nAv), so 


that 
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Averse 

between trackx 
Au. / Av 



Fig. 7 j Points on a rectangular grid in the (u,v) plan# at 
which the visibility is sampled for use with the 
discrete Fourier transform. A© shown, the spacing® Au 
and Av are equal. The division of the plane into cells 
for interpolation using the c el 1 -averaging technique 
is also shown. 
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Fig* 2*0 : Rearrangement of a visibility function for input to 
the discrete Fourier transform. The (u,v) origin Is at 
0, The contours represent the unsampled visibility 
which is measured in the range of u and v within the 
area ABCD. In the transformation of the sampled data 
the function is, in effect, replicated at points such 
as P, Q, and R. The input and output to the transfor- 
mation process cover only the positive values of u and 
v and are represented by the area OPQR. 
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£mAu , nAvj 


oo 




B(kA^ , l Ay? ) 


~j2n(kmAuA£ + LnAvAy>) 
e A£A n 


and 


B 



oo 



00 

V(mAu , nAv ) 

n=-*oo 



j 2rc ( kmAuA£ + 1 nAv Ay; ) 
e AuAv 


The Fourier transform relationship between the sampled 
functions requires that they be regarded as periodic, with 
periods MAu and NAv in u and v and NA£ and NAy> in £ and y> . This 
is illustrated in Fig. 2.8. 


A very important problem that has to be addressed is the 
problem of obtaining values of V(u,v) at the required points in 
the (u , v) plane. These points form a rectangular grid, as shown 
in Fig. 2.7, while the measured visibility values lie on the 
elliptical spacing locii. A very simple means of obtaining values 
at the grid points is to divide the (u,v) plane into rectangular 
cells of dimensions Au by Av centred on the grid points as shown 
in Fig. 2.7. The mean of all measured points within each cell is 
then computed and assigned to the grid point at the cell center. 
Grid points for which there are no measurements within the cell 
are assigned a value of zero. This procedure has been termed cell 
averaging. A 15 dB gaussian taper is applied to the visibility 
data in the grid. This reduces the aliasing of false details from 
outside the region for which the mapping is being done. 
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2.6 MUTUAL COHERENCE FUNCTION 

The generalised mutual coherence function is the measure of 

correlation in a radiation field between two points r. and r at 

1 2 

a time difference r and is defined as 

r (r l’ r 2 ’ a ^ = < S ( r 1 > S *(* 2 ’ t+T > > 

where S(r,t) is the response of the probe which is proportional 
to the radiation field E( r , t ) at any point in space r at any time 
t. For the special case of t=o 

r ( *V ^2 ’ 0) = r ( ^i- ‘V = < S(r lP t) S*(r 2 ,t)> 

is known as the spatial coherence function of the field. This 
coherence function when measured with the limitations of the 
finite bandwidth and beamwidth of the amplifiers and the antenna 
is known as the Visibility function, defined earlier. As has been 
mentioned earlier the relationship between the spatial coherence 
function and the brightness distribution is given by the Van 
Citterte Zernicke theorem, a simplified version of which is 

T(u, v) BC? . r/) 

If r be the discretized and sampled MCF and T be the discretized 

m 

and sampled MCF with the effect of the ionosphere then we can 
define the following Fourier relationship 
T(u,v) U(u , v) 44 B({ , r/) 

r m (u,v) U(u,v) 44 B^C? , rj) 


and 
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where U(u,v) is the sampling function which depends upon the grid 
points sampled in the (u-v) plane as explained earlier. 

Also 

V'W “ r(u i’ v i ) I(u i-V 

where I ( u^ , ) defines the effect of the ionosphere on the 
measured spatial coherence function. 

Henceforth in the thesis the terms visibility function and 
the spatial coherence function have been used interchangeably. 
The factor I(u,v) will be shown to depend upon the structure 
function of phase fluctuations in Chapter 4. 



CHAPTER - 3 


IONOSPHERIC SCINTILLATION AND MUTUAL COHERENCE FUNCTION 

Aa described in Chapter 1, the variation in the amplitude 
and phase of the signal is known as scintillation of amplitude 
and phase respectively. Most of what follows in the next 3 
sections is based upon the analysis presented in (Yeh and Liu, 
1982), (Ishimaru, 1978) and (Tatarskii, 1971). 

3.1 IONOSPHERIC IRREGULARITIES AND THEIR CHARACTERIZATION 

As has been noted earlier in Chapter 1, the permittivity of 
the ionosphere is given as 


. 2 2 
4n v £■ m 
o 

where N -> electron concentration 
e 

e -> charge on an electron 
and m -> mass of the electron. 

Due to random variations in the density of electron in the 
ionosphere, the emerging wavefront experiences a change of 
optical phase given by 

A<jfr = J An dl 

path through 
ionosphere 

where An is the refractive index fluctuation caused by electron 
density fluctuatios An. This change in the wavefront phase is 
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Fig. 3. 1 : The geometry of the problem : The thin phase screen 

approximation. The position vector of a point is 
defined as r = x i ^ + yi^ + Zl z = P + ZL Z 
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what causes scintillation, both phase and amplitude. To study 
further the effects of the ionosphere we consider the thin screen 
model of ionospheric irregularities. In this model, the 
irregularities are confined to a screen extending from z=o to 
z = L. The geometry is given in Fig. 3.1. Now we shall define some 
quantities of interest to be used later. 

1. Correlation function and spectra 

If AN(r) is the fluctuation of electron density from the 
background N q , then if AN ( r ) is assumed to be a homogeneous 
random field with a zero mean, then its correlation function Is 
( r ^ - r £ ) = <AN(r^) AN ( r^ ) > , where the angular brackets denote 
the process of ensemble averaging. By the Ui ener-Khinchin 

theorem, the correlation and the spectrum form a Fourier 
transform pair. 

oo - - 

j K i? 

*AN CK) = (2„r 3 JJ B an (F) . ' d 3 r 

~06 

B AN (?5 - JJJ « J " d3K 

If the irregularities are Isotropic, the correlation function 
depends only on Jr^-r^J. In-situ measurements and scintillation 
studies have revealed a variation of the form 

* AN (K) a K _p 


where isotropic! ty of irregularities is assumed. 
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The values of p have been found to be close to 4. However, 
in reality it has been found that the irregularities are 
anisotropic and aligned along the earth’s magnetic field. In the 
work of Rino and Fremouw (1977) it is shown that for a fairly 
general anisotropic model 

*AN Cq) = ab < An2 >-QCO 

where q 2 = AK 2 + BK K + CK 2 
x x y y 

The parameters a and b are axial ratios along and transverse to 
the principal irregularity axis, respectively. 

The function Q ( q) gives the shape of the spectral density 
function (SDF). It is normalized to unit area. Thus for a power 
law SDF, the appropriate functional form is 

8 n 3/2 r(v+l/2 )/r(r*-l) q 2v " 2 
Q(q) = ° 

where is the outer-scale cut-off wave number. The coefficients 
A, B and C depend on the propagation angles relative to the 
principle irregularity axis, and r (u- ) is the gamma function. 

The form of SDF is given in Fig. 3.2. By making the 
def ini t ion 

C = 8 n 3/2 < AN 2 > q 2v ' 2 T(v+l/2) / r(t>-l) 
s o 

it follows that 

abC ”... 

4 AN Cq:) = T 2.V+1/2 ' 

(q o +q ) 





Fig. 3.2 i The spectral density function Q(q) 
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3.2 OPTICAL PATH STRUCTURE FUNCTION 

At times the electron, density fluctuat ions end hence the 
optical path fluctuation as defined by 

<50 ( p) = J An (p, z) dz 

where p = (x,y) is the transverse co-ordinate and An is the 
refractive index fluctuations caused by electron density 
fluctuations . 


There is a background trend as they are not strictly 

homogenous but only locally homogenous. In these cases it is more 

convenient to deal with the structure function D defined by 

[ 2 
AN ( r + r’) - AN(r ’ )1 > 

The structure function for the optical path D, .ip') is just the 

A <p 

mean square value of the optical path difference between two 
points separated by p on the z = constant plane. This optical 
path structure function can be shown to be 




<|<50 (p+p’) - <50 (p')J > 

2 £< < 5 2 0 (p’) > - < <50 (p’ ) <50 (p + p') >J 

2 [% <°> - % r ^] 


where R^( - p 2 ) = < 60(p 1 ) <50 (p 2 ) > 
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The form R has been given in (Rino, 1979 ) and the result is 



(y) 


2 2 

r X L sec & G C 
e s 


, y ,v-l/2 

* 2q • 
o 


2»r(w + I) 


which gives 


D^(Ap) = 


2 2 

r X (C L) sec & G 
— % 1/2 '*' • '•(*’ - 


2 nr O + ±) 


2K j>-l/2 Cq o y) C 2 


y ,v-l/2 ] 

i > J 


where 


X 

C 

i 

q. 


-* classical electron radius 

-+ wavelength of incident radiation 

•* strength of turbulence as defined earlier 

-* spatial frequency corresponding to outer scale size 


K^(x) -+ modified bessel function of second kind and order v 
G -+ A geometrical factor 


C A p 


sx 


„ . . , 2 , 1/2 

BA p p + A A p 

sx sy sy 


AC 


B 

T~ 


where A, B and C are factors depending upon the angle of 

incidence and anisotropy and Sp and A p are quantities related 

sx sy 

to the differences in coordinates of the observing point, Ap. 

3.3 MCF AND THE STRUCTURE FUNCTION 


In this section the r elat ion between MCF and the structure 
function is explained. For this we consider the simplified case 



of normal incidence. The geometry of propagation is the same as 
given in chapter 2. 

In this case a component of the wave electric field is 
usually expressed as 

E(r,t) = u(r,t) exp £-j (kz - vt)j (3.1) 

where u(r,t) is the complex amplitude of the wave field and k and 
w are the wave number and the angular frequency of the wave. The 
time scale for temporal variation of the irregularities is 

such that w >> 2n /r . Hence under the forward scattering 
assumption the complex amplitude u still satisfies the parabolic 
equation 

" 2jk WZ + = -k 2 ^ (r f t )u (3.2) 

where represents the fluctuating part of the refractive index 
of the medium relative to the average refractive index. 


It has been demonstrated (Lee, 1974) that the MCF, T 
satisfies the following equation 



- (4nr 2 /k 2 ) |^ AN (0,0) - A Apj (p - p' , t-t ’ )J V (3.3) 

where is the classical electron radius and * T ) the 
integrated space-time correlation function for the electron 
density fluctuations AN(r,t) : 
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00 

A AN C *' t) = J B aN z > T) dz (3.4) 

~0D 

On making the co~ordinate transformation ( = p—p^ f X = (p+p ’ ) / 2 
and using the fact that r must be independent of X for a 
homogenous random medium, it is seen that r(jf, z) satisfies the 
equations 


- mL = 

2, 2 

r l 


dz 

-r X 
e 

&r 

[ A AN C0,0) " A AN a ’ T >| r < OSz^L 

(3.5) 

and 

= o 

dz 

; z > L 



For the spatial coherence function r=0 and the solution for the 
above equations becomes 


TCAp.z) = r(Ap.O) exp -r ^ \ 2 L |A an ( 0,0) - A^CAp 


which can be rewritten as 


H A [- 


-■>]) 


r(Ap,z) = r (Ap,0) exp 


- v CAp) 


(3.6) 


or 


T(u,v) = r (u,v) exp 


°A* 


(3.7) 


r ~D (ApX 

The term exp J — | 1 defines the ionospheric contribution 

term I(u,v). Even though the above relation has been derived for 
normal incidence, it is true for generalized oblique incidence 
also (Rlno, 1978). 



CHAPTER - 4- 


PROGRAM CODE DESCRIPTION 

The aim of this work is to study the effects of the 
Ionosphere on reconstructed radio-images. The program code can be 
seperated Into two distinct parts. In fact the first part 
makes use of the second. 

The main program is the first and this part is responsible 

for 

i) Calculating the U(u,v) function 

ii) Generating the model gaussian sources 

iii) Calculating the various forward and inverse transforms, 
and iv) Calculating the maximum and r.m.s. dB errors between 

the assumed source image (A) , the reconstructed image 
without the effect of ionosphere (B) and the 

reconstructed image with the effect of the ionosphere 
(C) . 

U(u,v) function behaves as a spatial frequency filter. It 

arises due to the inherent limitations of the process of 

measurement of spatial frequencies by the synthesis telescope. 
This is introduced due to the finite size of the telescopic 

array. If we could somehow measure all the (u,v) frequency 
components then this would be a constant function. Since this is 
not so, we have a two-dimensional function with ’holes’ in it 
corresponding to the unsampled frequency components. These 
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b l 

= (D/X). 

( y i ~ y 2 ) 

. cos 

4 > 





(4. 

1) 


b 2 

- (D/X). 

(y rV 

. sin 

4 > 





(4. 

2) 


u 

= b^ sin 

H 







(4. 

3) 


V 

= b^ cos 

6 - b 2 

sin 6 

cos 

H 




(4. 

4) 


where S = declination of the source, 

H = hour angle of the source, 

4> - latitude of the telescopic site 
D = baseline length 


(x,y,z) define the coordinates of antennas forming the 
interferometer, with (x,y) being measured in a plane parallel to 
the earth’s equator, x in the meridian plane, y towards the east 
and z is measured towards the celestial north pole (Fig. 2.6). 
The tracking of source is done for all the time the source has an 
elevation angle greater than fifteen degrees. The calculations 
for u and v are done after intervals of 1 min. This interval has 
been fixed by taking into account the speed of coverage of the 
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u-v plane by the u-v ellipses. A compromise has to be worked 
between the facts that more frequent samples may result in the 
points lying in the same grid whereas less frequent calculation 
may end up missing some intermediate grid points through which 
the ellipse passes. A 128 x 256 array is used to store the number 
of times each (u,v) cell is sampled. The cells which are not 
sampled are assigned a value of zero. Use is made of the fact 
that , 

V( u , v ) = V* ( -u , -v) 

So whenever a point ( u , v ) is sampled, the corresponding point 
(-u, -v) is also sampled. 

The major source of distortion in reconstructed Images is 
the existence of unsampled cells in the u-v plane. The U(u,v) 
function tells us the extent of u-v plane coverage. Fig. 4.1 
gives us an Idea of the U(u,v) function for different 
declinations and frequencies for the GMRT 30 element array. The 
configuration of the GMRT array is given in Fig. 4.2 (Swarup, 
1990) and the co-ordinates of the various antennas are given in 
Table 4.1. 

The second part is a routine for the calculation of the 
structure function which is used for calculation of the 
ionospheric effect I(u,v) as explained earlier. 









TABLE 4.1 


The GMRT antenna array co-ordinates in metres 


Arxt enna 

X 

i . 

-258.0 

2 . 

802.0 

3 . 

103.0 

4 . 

237.0 

5 . 

998.0 

6 . 

360 . 0 

7 . 

461 . 0 

8 . 

149.0 

9 . 

386 . 0 

10. 

600.0 

11. 

1034 . 0 

12. 

266 . 0 

13. 

847 . 0 

14. 

903.0 

15. 

-2690.0 

16. 

11610.0 

17 . 

-4190.0 

18 . 

-7590.0 

19. 

-10090.0 

20. 

-11890 . 0 

21 . 

-290.0 

22 . 

760.0 

23. 

-240.0 

24 . 

-790.0 

25 . 

660.0 

26 . 

1810.0 

27 . 

3510.0 

28 . 

5710.0 

29 . 

7310.0 

30 . 

8410 . 0 


Y 

1 z 

259 . 0 

33.0 

94.0 

39.0 

237 . 0 

35.0 

259.0 

37.0 

104.0 

40.0 

497.0 

42.0 

469.0 

43.0 

659 . 0 

49 . 0 

396 . 0 

44.0 

850.0 

45.0 

575.0 

47.0 

907.0 

15.0 

1380.0 

17.0 

899.0 

19.0 

- 800.0 

53.0 

- 9250.0 

01 . 0 

- 1750.0 

55.0 

- 2750 . 0 

58.0 

- 3350 . 0 

61.0 

- 4300 . 0 

63.0 

2750 . 0 

31.0 

4950.0 

21.0 

7150.0 

29.0 

9750.0 

23.0 

14400.0 

25.0 

- 300.0 

13.0 

- 1300.0 

11.0 

- 2600.0 

09.0 

- 5200.0 

07 . 0 

- 8000.0 

05.0 
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The aubroutine ION calculates the factor 
I(u,v) = exp ( -D</>/ 2 ) 

This subroutine makes use of the anisotropic ionospheric 
irregularity model developed by Fremouv and Secan (1984) for 
their UBMOD program for calculation of ionospheric scintillation 
indices. The inputs to the program required are : 

i) Time of the day (Diurnal effect) 

ii) Elevation angle 

lii) Azimuth angle 

iv) Sun spot number 

v) The planetary magnetic index Kp. 

For relating the time of the day to the hour angle of the source 
we have made use of the relations 

Local civil time (LCT) = Local sidereal time (LST) - x 4 M x D 

(4.5) 

LST = R.A. of source (ot) + Local hour angle (LHA) (4.6) 

so that 

LCT = a + LHA - ir x ^ D (hrs.) (4.7) 

"t hi 

where D = Day of the year reakoned from the 264.7 day of the 

t 

year (September 21 ) which is the autumnal equinox. 
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Using the standard relations for relating the azimuth angle 
and the elevation angle to the hour angle, we can for each sample 
calculate the value of I(u,v). A 256 x 256 array (KDATA) stores 
the value of I(u,v) in the proper grid position. For each new 
sample corresponding to the (u.^ , v ) point, the value in that 
KDATA cell is incremented by I( u 1 , v 1 ). After the completion of 
the tracking, the average value of the measured visibility 
function is calculated by dividing the values of KDATA array by 
the corresponding total number of samples. The unsampled points 
are assigned a value of zero. 

Uhen there is no ionospheric effect, we have a 256 x 256 
array filled with l’s and 0's corresponding to the sampled and 
unsampled points. Now a uniform gaussian taper of 15 dB at the 
edges is applied and the inverse FFT is calculated to get the 
reconstructed image from the Fourier transform of the original 
assumed source distribution. 

For the KDATA array, the same procedure leaves us with an 
array filled with zeroes corresponding to unsampled points and 
values different from zero for the sampled points. After applying 
the gaussian taper and performing the inverse FFT we get the 
reconstructed image with the effect of the ionosphere. 

The dB errors - maximum and r.m.s. for Images A, B and C are 
calculated for different regions surrounding the source. For the 
purpose of calculations a point source and an extended source 
have been considered. 



CHAPTER 5 


RESULTS AND CONCLUSION 

The effect of following parameters on the reconstructed 
image were studied 

a) The size of the assumed source brightness distribution 

b) Declination of the source 

c) Frequency of observation 

d) Right ascension of source 

e) Solar activity in terms of sun spot number 

f) The season of the year. 

Table 0 lists the figure and quantitative data table numbers 
for the different parameters for which study has been conducted. 
Note that Plot-A corresponds to the assumed source distribution, 
B to the reconstructed source distribution without the effect of 
ionosphere and C to the reconstructed source distribution with 
the effect of ionosphere. The Fig. 5.1 gives the Plot-A for the 
point source and Fig. 5.2 for an extended gaussian source with 
the 3 dB half width of 60 arc-seconds. The planetary magnetic 
index has been kept at a typical value of 4.0 since dependence on 
this parameter is very weak. 

5.1 EFFECT OF SIZE OF SOURCE 

Two types of source distributions were studied 

a) Point source brightness distribution 

b) Extended gaussian source brightness distribution. 



Table O : List of various figures and tables for different 
parameters 


Source 

SSN 

Point 


Extended 


Point 

50 

Point 

50 

Point 

50 

Point 

50 

Point 

50 

Point 

50 

Point 

50 

Point 

200 

Ext ended 

50 

Extended 

50 

Ext ended 

200 

Point 

50 

Point 

50 

Ext ended 

50 

Extended 

50 

Ext ended 

50 


Day <5 Freq. 

(degrees) (MHz ) | (hours ) 


Sun Spot Number 
Planetary Magnetic Index 

70 (March) 

170 (June) 

Right Ascension 

Declination 
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In general , spurious structure are introduced into the 
reconstructed images even without the effect of the ionosphere 
being considered. This spurious structure is introduced due to 
the truncation of the sample space based on maximum baseline 
available, measurement at discrete intervals and existence of 
unsampled points in the u-v plane. 

There is stretching of the reconstructed source image, 
roughly along constant right ascension axes. The reconstructed 
image which includes the effect of the ionosphere is broadened in 
a direction perpendicular to the direction of stretching of the 
reconstructed i stage without the effect of the ionosphere. This is 
borne out quite clearly in the Figs. 5.4 and 5.6. 

The rms distortion In an area of radius 30 arc seconds is 

less than 3 dBs for the point source and is less than 2 dB for 

the extended source even within an area of 60 arc-seconds. For an 
extended source, the distortion levels are very small. The 
distortion level as 0 dBs upto a radius of ^ 60 arc-seconds which 
is to be expected for a source with a 3 dB half width of — 60 

arc-seconds. For very severe cases of distortion one cannot even 

distinguish between a point source and extended source as is seen 
from Figs. 5.5 and 5.11 and Figs. 5.7 and 5.12. 
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5.2 EFFECT OF DECLINATION OF SOURCE 

Two effects are observed. 

i) The increased tracking time for higer declinations leads Lo 
greater fidelity reconstructed images 

ii) Tracking at large zenith angles for larger portion of Lime 
leads to greater distortion. 

These two effects are opposing in nature and are seen 
clearly from plots in Figs. 5.3, 5.4, 5.5 and 5.6, which 
correspond to declinations of -30, 0, 30 and 60, respectively. 

The better fidelity can be seen from smaller r.m.s. error 
between plots (A-B) for increasing declination and the greater- 
distortion can be seen from the increasing r.m.s. error between 
plots (B-C) and (A-C) as seen from Tables 1, 2, 3 and 4. 

5.3 EFFECT OF FREQUENCY 

As is to be expected, there is increasing distortion with 
decreasing frequency. The effect of ionosphere at a declination 
of 30° can be clearly seen in Figs. 5.5, 5.8, 5.14, 5.15 
(corresponding to frequencies of 38, 153, 232 and 327 MHz 
respectively) with Figs. 5.15 and 5.16 showing the comparative 
effects on point and extended source, respectively at 327 MHz . 
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5.-4 EFFECT OF RIGHT ASCENSION 

The RA determines the time of the day, during which the 
source is being tracked. Since the maximum scintillation occurs 
around 22:30 hours local time, we expect maximum distortion to 
take place for images which are produced by tracking around this 
time. This is infact true as is borne out by Figs. 5.5 and 5.9 
for a point source. Tables 3 and 7 also bring out these factors 
quantitatively . 

5.5 EFFECT OF SEASON 


The seasonal scintillation maxima occur during March 
and October and minima around December and June. Hence, we expect 
the results to change accordingly. Figs. 5.11, 5.17 and 5.18 
show the effect of the effect of season on an extended source at 
38 MHz. Figure 5.18 shows the effect of the season for the same 
parameters except inclusion of the night time maxima. There is 
greater stretching of 6 dB lobes along the constant RA axes. This 
implies greater distortion of the image, even though this may not 
be reflected in the corresponding quantitative figures. 


S. e EFFECT OF SUN SPOT NUMBER CSSNJ 


The magnetic activity is parameterised in terms of the sun 
spot number. There is no appreciable difference as is shown in 
Figs. 5.5 & 5.10 (for SSN = 50 and 200 at 153 MHz) and Figs. 5.11 
i 5.13 (for SSN - 50 and 200 at 38 MHz) together with the 
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correspond! aft quantitative errors in Tables 3, 8, 9 and 11. This 
is to be expected as there is very weak relationship between 
scintillation strength and magnetic activity. 

5. 7 CONCLUSI ON 

There would be large distortion for 38 MHz and 153 MHz 
reconstructed images. The images of sources, both point and 
extended will be extended, more or less in a direction parallel 
to axes of constant right ascension. For month of March and 
October, It will be very difficult to distinguish between point 
and extended source and we might as well discard data for higher- 
declinations for these months. This is true even for data which 
may not include the night-time scintillation maxima. 

There is only a small effect of scintillation on the 233 MHz 
and other higher frequency observations and we can safely neglect 
these effects. 





Figure 5. lCbl 


Fig. 5.^ * C aD Assumed point source distribution 
Cbl A blown out view of source in Cal 



Figure 5. a) 



Figure 5. 2CbI> 


Fig. 5. 2 i C al> Assumed extended gauss! an source distribution 
C b J A blown out view of source in C a! 
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Fig. 5. 3C aJ : Plot-B corresponding to Plot-A of' Fig 


Fig. 5. 3( b) s Plot —C corresponding to Plot —A of* Fig. 5. 1 


max 


Table 1 : dB Errors for Plots of Fig. 5.3 


A- 

-B 

max 

rms 

9.0 

3 . 7 

19 . 0 

10.9 

33.0 

20.9 

75.0 

47.7 

87.0 

50.4 
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Fig. 6. SC sO i Plot— B cor responding to Plot-A of Fig. 5.1. 



Fig, 5, 5C b) s Plot — C corresponding to Plot, —A of - Fig. 5, 1 « 


— 
— 1 

A- 

-B 

A- 

-c 

B- 

-C 

max 

rms 

max 

rms 

max 

rms 

0 

1 

u> 

o 

8.0 

3 . 6 

m 

3.8 

1.0 

0.4 

0 

1 

Os 

o 

18.0 

10.8 

bB9 

11 . 2 

2.0 

0.8 

0-120" 

32.0 

20.4 

199 

21.1 

2.0 

1.2 

0-300" 

72 . 0 

46.0 

ill 

48 . 6 

16.0 

4.5 

0-600” 

83.0 

48.3 

m 

53 . 3 

22.0 

9.3 


Table 3 


dB Errors for Plots of Fig. 5.5 
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t r ack i ng f a r l £ i £S h r s 
■f roroi *6: ?8hrs to: 4 : lOhrs 


'tracking far: i2*gg hr* 
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'$5N: m.m HP : ' 4.80' Dfty'i 170.® 
Dec; 60*90 Frq; 38.43 Rfl : 16.00 


' -r V- 


r-^8 «»6® «*■»* -ij® 


rms 


max 


rms 


max 


Table 4 


8.0 

3.5 

1.0 

17.0 

10.7 

1.0 

28.0 

19 . 9 

5.0 

75 . 0 

46.3 

16.0 

89 . 0 

54.5 

25.0 


5. 1. 
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Plot - ro 
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Fig. 5. 7Ca) 1 Plot-B corresponding to Plot-A of Fig 


Fig. 5. 7Ch'J i Plot — C corresponding to Plot-A of Fig 


rms 


max 


Table 5 : dB Errors for Plots of Fig. 5.7 


5.1. 
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Table 6 : dB Errors for Plots of Fig. 5.8 
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A- 

-B 

max 

rms 

0-30” 

8.0 

3.6 

0-60” 

18.0 

10.8 

0-120” 

32.0 

20 . 4 

0-300” 

72.0 

46.3 

0-600" 

83.0 

48.3 


B- 

-c 

max 

rms 























Fig. 5. UCal) : Plot-B corresponding to Plot-A of Fig. 5.2 



Fig. 5. 1 1 C b) : Plot-C corresponding to Plot-A of 


Fig. 5. &. 


Region 

A - 

-B 

A- 

-c 

B- 

•c 

max 

rms 

max 

rms 

max 

rms 

0-30" 

n 

0.0 

0.0 

0 . 0 

0 .'o 

0.0 

0-60” 

warn 

0 . 2 

1.0 

0 . 5 

1.0 

0 . 5 

0-120” 

1 . 0 

0.2 

1 . 0 

0 . 7 

1 . 0 

0.7 

0-300” 

1 . 0 

0 . 4 

3.0 

1 . 2 

3 . 0 

1.1 

0-600” 

1 . 0 

0 . 5 

5.0 

2 . 2 

5.0 

2 . 1 


Table Q : dB Errors for Plots of Fig. 5.11 
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Fig. 5 . USCa} : Plot-B correspond! ng to Plot -A of Fig 



Table 1 0 > dB Errors for Plots of Fig 


5 . 1 2 CbD : Plot-C corresponding to Plot -A of Fig. 
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g. 5. ISCb'J: Plot-C corresponding to Plot -A of Fig- 5-2 
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Fig. 5. 15 CjO i Plot-B corresponding to Plot-A of Fig 


Fig. 5. 15Cb!> t Plot— C corresponding to Plot-A of Fig 
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A- 

-B 

A- 

■C 

B- 

■C 
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max 
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5.0 
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5.1. 


Table 13 


dB Errors for Plots of Fig. 5.15 
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Fig. 5. 18Ca3 i Plot-B correspond! ng to Plot -A of Fig. 5.2. 



Fig. 5. 1 8C fcO s Plot-C corresponding to Plot-A of Fig. 5.2. 
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A- 
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A- 

-C 

B- 
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rms 
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0-30” 
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1 . 0 
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0-120" 

1 . 0 


1.0 


1 . 0 

0 . 6 

0-300” 

1.0 

0 . 4 

3 .0 


3 . 0 

1 . 1 

0-600" 

1 . 0 

0 . 5 

5 . 0 


5 . 0 

2.0 


Table 1 6 i dB Errors for Plots of Fig. 5.18 
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